Quark-Hadron Duality in Neutron ( 3 He) Spin Structure 

P. Solvignon 21,1 , N. Liyanage 22 , J.-P. Chen 12 , Seonho Choi 21 , K. Aniol 3 , T. Averett 23 , 
W. Boeglin 7 , A. Camsonne 2 , CD. Cates 22 , C.C. Chang 15 , E. Chudakov 12 , B. Craver 22 , 
F. Cusanno 9 , A. Deur 22 , D. Dutta 4 , R. Ent 12 , R. Feuerbach 12 , S. Frullani 9 , H. Gao 4 , 
F. Garibaldi 9 , R. Gilman 12,19 , C. Glashausser 19 , V. Gorbenko 11 , O. Hansen 12 , 
D.W. Higinbotham 12 , H. Ibrahim 18 , X. Jiang 19 , M. Jones 12 , A. Kelleher 23 , 
J. Kelly 15 , C. Keppel 8 ' 12 , W. Kim 14 , W. Korsch 13 , K. Kramer 23 , G. Kumbartzki 18 , 
J.J. LeRose 12 , R. Lindgren 22 , B. Ma 17 , D.J. Margaziotis 3 , P. Markowitz 7 , 
K. McCormick 18 , Z.-E. Meziani 20 , R. Michaels 12 , B. Moffit 23 , P. Monaghan 17 , 

C. Munoz Camacho 5 , K. Paschke 16 , B. Reitz 12 , A. Saha 12 , R. Sheyor 20 , 
J. Singh 22 , K. Slifer 21 , V. Sulkosky 23 , A. Tobias 22 , CM. Urciuoli 10 , K. Wang 22 , 
K. Wijesooriya 4 , B. Wojtsekhowski 12 , S. Woo 14 , J.-C. Yang 6 , X. Zheng 1 , L. Zhu 17 

(Jefferson Lab E01-012 Collaboration) 

1 Argonne National Laboratory, Argonne, IL 60439 
2 Umversite Blaise Pascal et CNRS/IN2P3 LPC, 63177 Aubiere Cedex, France 
3 California State University, Los Angeles, Los Angeles, CA 90032 
A Duke University, Durham, NC 27708 
5 CEA Saclay, DAPNIA/SPhN, F -91191 Gif sur Yvette, France 
6 Chungnam National University, Taejon 305-764, Korea 
7 Florida International University, Miami, FL 33199 
8 Hampton University, Hampton, Virginia 23187 
9 Istituto Nazionale di Fisica Nucleare, 
Gruppo Collegato Sanitd, Seziona di Roma, 00161 Roma, Italy 
10 Istituto Nazionale di Fisica Nucleare, 
Sezione di Roma, 00185 Roma, Italy 
n Kharkov Institute of Physics and Technology, Kharkov 61108, Ukraine 
12 Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 
13 University of Kentucky, Lexington, KY 40506 



1 



14 Kyungpook National University, Taegu City, South Korea 

15 University of Maryland, College Park, MD 20742 

16 University of Massachusetts, Amherst, MA 01003 

17 Massachusetts Institute of Technology, Cambridge, MA 02139 
18 Old Dominion University, Norfolk, VA 23529 
19 Rutgers, The State University of New Jersey, Piscataway, NJ 08855 
20 University of Tel Aviv, Tel Aviv, 69978 Israel 
21 Temple University, Philadelphia, PA 19122 
22 University of Virginia, Charlottesville, VA 22904 
23 College of William and Mary, Williamsburg, VA 23187 
(Dated: October 26, 2008) 

Abstract 

We present experimental results of the first high-precision test of quark-hadron duality in the 
spin-structure function g\ of the neutron and 3 He using a polarized 3 He target in the four- 
momentum-transfer-squared range from 0.7 to 4.0 (GeV/c) 2 . Global duality is observed for the 
spin-structure function g\ down to at least Q 2 = 1.8 (GeV/c) 2 in both targets. We have also 
formed the photon-nucleon asymmetry A\ in the resonance region for 3 He and found no strong 
Q 2 -dependence above 2.2 (GeV/c) 2 . 

PACS numbers: 13.60.Hb, 13.88,+c, 14.20.Dh 
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Quark-hadron duality describes the remarkable similarity between electron-nucleon scat- 
tering in the deep inelastic scattering (DIS) region, where the electron scatters off an asymp- 
totically free point-like quark, and in the nucleon resonance region where the electron scatters 
off a highly correlated cluster of quarks and gluons. 

In high-energy electron-nucleon scattering, the electron interacts with the nucleon con- 
stituents by exchanging a virtual photon. This process depends on two variables, Q 2 and 
x. The quantity Q 2 is the four-momentum squared of the exchanged virtual photon and 
x = Q 2 /(2Mu) is the Bjorken scaling variable, which represents the fraction of the nucleon 
four-momentum carried by the struck parton, with M the mass of the nucleon and v the 
energy transfer between the lepton and the nucleon. The invariant mass W of the virtual 
photon- nucleon system is related to Q 2 and x through W 2 = M 2 + Q 2 /x — Q 2 . The scattering 
cross section is parametrized in terms of two unpolarized structure functions (iq and F%) and 
two polarized structure functions (gi and g?). These structure functions are related to the 
polarized and unpolarized quark distribution functions in the nucleon. Another observable 
that provides a direct insight into the polarized quark distributions in the nucleon is the 
photon-nucleon asymmetry, A±, which is approximately equal to the ratio of the polarized 
and unpolarized structure functions, pi/iq, at large Q 2 . 

In the DIS region, at asymptotically large values of Q 2 and W, the electron scatters 
off a point-like free quark. As a result, the structure functions are independent of Q 2 (en- 
sealing), modulo logarithmic corrections due to gluon radiation. On the other hand, at low 
values of Q 2 and W, the quarks strongly interact with each other and resonance production 
dominates. In this nucleon resonance region, the electron can be considered to be scattering 
off a correlated cluster of quarks and gluons, the nucleon. Thus, the structure functions are 
expected to depend strongly on Q 2 . However, data taken at Stanford Linear Accelerator 
Center (SLAC) more than 30 years ago revealed that the scaling curve seen for DIS data 
at high Q 2 is an average over the resonances at lower Q 2 but at the same value of x 
Since then, there have been many experiments to study this intriguing observation and 
numerous theoretical attempts to explain it. For a review of the full body of experimental 
and theoretical work see Ref. 

In the case of unpolarized electron scattering, duality can be expressed as 

X2(W2,Q 2 ) I-X2 

dx F 2 res (x,Q 2 ) = / dx F 2 DIS (x,Q 2 ), (1) 



xx(Wx,Q 2 ) J x x 
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where _F| es (;c, Q 2 ) is the structure function measured in the resonance region at low Q 2 while 
F^ ls (x,Q 2 ) is the scaling function evolved down to the same Q 2 . Global duality is said 
to hold when Eq. [T]) is performed over the entire resonance region. Local duality would be 
observed if the above equality holds when integrated over a resonance bump jj, 4]. Data from 
Jefferson Lab (JLab) Hall C [5j have demonstrated that for unpolarized structure functions 
Bloom-Gilman duality holds at the 10% level down to Q 2 = 0.5 (GeV/c) 2 . These data also 
show that local duality holds for each of the three prominent resonance regions. 

The observation of duality in unpolarized structure functions raises the exciting question 
whether quark-hadron duality holds for spin-structure functions as well. Indeed, duality 
in spin-structure functions would be even more intriguing than in the unpolarized case, 
since the spin-structure functions come from the differences of cross sections, they are not 
necessarily positive definite. As a result, a spin-structure function measured at the same 
value of x in resonance and DIS regions can, in principle, have opposite signs. For example, 
the contribution to g\ at the A resonance is negative at low Q 2 , while the scaling curve for 
g\ is positive. 



Recent data from Jefferson Lab 



2j and DESY [8J have shown that quark-hadron duality 
holds globally, in the case of polarized structure functions of the proton and the deuteron, 
down to Q 2 = 1.7 (GeV/c) 2 . However, the high-precision results from Jefferson Lab Hall 
B indicated that local duality is violated for the proton and deuteron polarized structure 
functions in the region of the A resonance even for Q 2 values as high as 5.0 (GeV/c) 2 js]. 

Polarized structure functions measured previously over the low Q 2 -range of 0.1-0.9 
(GeV/c) 2 , using a polarized 3 He target as an effective neutron target have shown hints 
of quark-hadron duality for the polarized structure functions of the neutron. However, until 
now there have been no resonance-region neutron spin-structure data in the intermediate 
Q 2 -range where quark-hadron duality is expected to manifest. 

In this letter we report results from Jefferson lab experiment E01-012 where we measured 
polarized structure functions in the resonance region over the Q 2 -range 0.7-4.0 (GeV/c) 2 
using a polarized 3 He target as an effective polarized neutron target. These new results 
provide a precision test of quark-hadron duality for neutron spin-structure functions. 

E01-012 took data in 2003 in Hall A at Jefferson Lab where inclusive scattering of longi- 
tudinally polarized electrons from a longitudinally or transversely polarized 3 He target was 
studied. We formed asymmetries and polarized cross section differences from data taken at 
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a scattering angle of 25° for three incident beam energies, 3.028, 4.018 and 5.009 GeV, and 
at 32° for an incident beam energy of 5.009 GeV. 

The polarized electron beam was produced from a strained GaAs photocathode which 
was illuminated by circularly polarized light. The helicity of the beam was pseudo-randomly 
reversed at a rate of 30 Hz. The beam polarization was monitored by a M0ller polarimeter 
and was found to be stable during each incident energy running but varied between 70 and 
85% for the production data, depending on the incident beam energy values and the status 
of the other experimental halls. The total relative uncertainty on the beam polarization was 
3.4%. The beam current delivered on the target was kept below 15 //A in order to limit 
target depolarization. 

The polarized target consists of two connected chambers with an admixture of rubidium 
(Rb) and 3 He in the upper chamber where spin-exchange optical pumping [10] takes place. 
About 90 W of 795 nm light produced by fiber-coupled laser diode arrays and an oven heated 
to 170°C were used to this goal. Under running conditions, the 3 He density was 3.20 x 10 20 
cm -3 in the lower chamber where the electron beam passed through and from which the Rb 
was absent due to the lower temperature. The polarization of the target was monitored by 
two independent polarimetries: one based on nuclear magnetic resonance (NMR) [ljj and the 



other on the electron paramagnetic resonance frequ ency shift (EPR) [12| . Both polarimetries 
rely on the adiabatic fast passage (AFP) technique ll]]. The average target polarization was 
0.38 ± 0.02. The cells also contain a small amount of nitrogen gas (at about 1% of the 3 He 
density) to quench radiative transitions which limit the optical pumping efficiency However, 
this introduces a dilution factor in the extracted physics quantities. Nitrogen dilution factors 
were determined from measurements on a reference cell filled with nitrogen gas and were 
found to dilute the asymmetries by 5 to 9%. The asymmetry of elastic scattering e — 3 He was 
extracted from data taken at an incident beam energy of 1.046 GeV and a scattering angle 
of 16° in order to obtain an independent measurement of the product of beam and target 
polarizations and nitrogen dilution. Theoretical and experimental elastic asymmetries were 
found to agree within a total uncertainty of 5% 13] . 



Both Hall A high resolution spectrometers 14| were set at the same momentum and angle 
configurations and tuned to optimize the detection of the scattered electrons. This symmetric 
configuration has the advantage of obtaining two independent measurements and therefore 
controlling part of our systematic uncertainties. A gas Cerenkov counter combined with a 
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two-layered lead-glass calorimeter provided particle identification with a pion rejection better 
than 10 3 to 1 while keeping the electron efficiency above 99%. The remaining contamination 
from pions to the asymmetries has a maximum effect on the asymmetries of 15 ppm I3I ]. 
which is negligible compared to the statistical uncertainties. 

We formed the longitudinal and transverse polarized cross section differences using: 

d 2 a^ d 2 a^ 

A(T|I(±) = -dndW - -dmEr (2) 

where the superscript j (j) represents the incident electron spin direction and ff the 
longitudinal (transverse) target spin direction. The quantities E' and Q correspond to the 
scattered electron energy and the spectrometer solid angle, respectively. The polarized 
structure function g\ can be formed directly from the polarized cross section differences as 
follows: 

* = 4af ^ + g ) |Ag,+tan 2 Agj (3) 

Here, a e is the fine structure constant, 9 is the scattering angle and E is the energy of the 
incident electron beam. 

External and internal radiative corrections have been applied to the polarized cross section 
differences following the formalism of Mo and Tsai 15| for the spin-independent part and 
the approach of Akushevich and Shumeiko [16] for the spin-dependent part. 

The structure function gi was generated at constant energies and scattering angles. Our 
goal in this analysis is the quantitative study of quark-hadron duality for the spin-structure 
function g\. Following the same procedure as in the unpolarized case (see Eq. [TJ, g\ would 
be integrated at constant Q 2 over a limited x-range. However forming this partial mo- 
ment requires the interpolation of g\ to a constant Q 2 . The coverage of our data allows 
interpolation of g\ to the four Q 2 -values of 1.2, 1.8, 2.4 and 3.0 (GeV/c) 2 . The results for 
the 3 He spin-structure function g\ are shown in Fig. [TJ The DIS parameterizations of q-\ 
for the proton and the neutron, from Bliimlein and Bottcher (BB) 



AAC collaboration [19( and Leader, Sidorov and Stamenov (LSS) 2G] 



17[, GRSV group (l8|, 



were evolved to those 



Q 2 -values and then combined to obtain the 3 He DIS parametrizations using the effective 
polarization equation [2l|: 

s! He = Vn9l + 2p p gl (4) 
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The effective polarizations of the neutron and the proton are p n = 0.86 ± 0.02 and p p = 
— 0.028±0.004, respectively j^]. While neglecting Fermi motion and nucleon binding effects 
(EMC effect), this method of computing the 3 He structure function from the proton and 
the neutron structure functions was confirmed to be a good approximation at the 4% level, 
which is well below the DIS parametrization precision. All parameterizations were taken at 
Next-to-Leading Order (NLO). Target-mass corrections [23| have been applied to the DIS 
parameterizations in order to take into account the finite mass of the nucleon. Also plotted 
are the highest precision world data on gl He in the DIS region from SLAC E154 2^] and JLab 
E99-117 25| . It can be observed that our resonance data at Q 2 = 1.2 (GeV/c) 2 oscillate 
around the DIS curves and are significantly lower in the A(1232) region (x ~ 0.65). As Q 2 
increases, our data approach the DIS parameterizations indicating quark-hadron duality. 

To perform a quantitative test of quark-hadron duality, g\ was integrated over the x- 
interval corresponding to the resonance region 26] covered by our data, which is from pion 
threshold (x w ) to an x corresponding to W = 1.905 GeV (x min ): 



fi(Q 2 



dx gi(x, 



(5) 



The experimental gi-integral for the neutron was extracted using the method described 
in Ref. [27J: 



Pn Pn 



(6) 



Data from JLab Hall B experiment on the proton spin-structure function 6| were used for 
rf. Figure [2] shows the comparison of the integral of g\ over the resonance region to the DIS 
parameterizations. Data from JLab E94-010 9J are also plotted and a smooth transition 
can be observed with our data. 

At the lowest Q 2 , our data deviate slightly from the DIS parameterizations for both 3 He 
and the neutron. As Q 2 increases, the resonance data and the DIS parametrization show 
good agreement indicating global quark-hadron duality holds at least down to Q 2 = 1.8 
(GeV/c) 2 and above. It is important to note that, for the proton and the deuteron, global 



duality was previously observed for Q 2 above 1.7 (GeV/c) 2 p 

Quark-hadron duality was also studied for the photon- nucleon asymmetry A\. From the 
parallel and perpendicular asymmetries A\ was extracted as follows: 

A {l (x,Q 2 ) r]A ± (x,Q 2 ^ 
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where 



1 + eR(x, Q 2 ) V 1 + e 2e 

_ e^/Q 2 = i 

77 E — E'e ' 1 + 2(1 + z/ 2 /Q 2 ) tan 2 § 

.D is the photon depolarization factor and e is the photon polarization. The longitudinal 
to transverse cross section ratio R(x, Q 2 ) has never been measured on 3 He in the resonance 
region. Precise data on this ratio exist for the proton and the deuteron in the DIS region 29 1 
and show a reduction of deuteron ratio by 30% compared to the proton, i.e. Rd = 0.7R P . 
Considering that quark-hadron duality holds for R the same reduction of the ratio R 
can be expected in the resonance region. Besides the ratio R is believed to diminish further 
for 3 He. Therefore, it was estimated for 3 He that R = 0.5R P with R p extracted from fits 
to the resonance data [30;]. We then studied the sensitivity of A\ to different models and 
assumptions for R. The systematic effects on A± were found to be significantly smaller than 
the relative statistical uncertainties and were added to the systematic uncertainties. 

The photon-nucleon asymmetry A\ for 3 He in the resonance region is presented in Fig. [31 
Also plotted are the DIS data in order to provide a direct comparison between the deep 
inelastic scattering behavior and the resonance data trend. We performed a parametric fit 
to the 3 He DIS data resulting in: x a580 (-0.135 + 0.336x - 0.174x 2 )(l + 0.578/Q 2 ). For 
Q 2 below 2.0 (GeV/c) 2 , it can be seen that A 3 ^ in the vicinity of the A(1232) peak is 
large and negative unlike the DIS behavior. At higher Q 2 , our A 3 ^ data cross zero and 
become positive even in the A(1232) region. This is due to the increasing importance of 
the non-resonant background with respect to the resonance strength. The most noticeable 
feature is that, for values of Q 2 above 2.2 (GeV/c) 2 , A^ He from our two data sets agrees 
very well indicating little or no Q 2 -dependence. Furthermore, the trend of these data show 
the same particularity as the DIS data from E99-117 of becoming positive at high x, a 



3JJ and Relativistic Constituent 



prediction for DIS data from both pQCD inspired models 
Quark models [32] ■ 

In summary, the results presented in this paper are the first precision test of quark- 
hadron duality for the neutron and 3 He spin-structure functions. Global duality has been 
demonstrated to hold for the neutron and 3 He polarized structure function g± down to at 
least Q 2 = 1.8 (GeV/c) 2 . We also observed no strong Q 2 -dependence in the photon-nucleon 
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asymmetry A 1 He measured in the resonance region for Q 2 > 2.2 (GeV/c) 2 and that A 1 He 
becomes positive at large x. 
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FIG. 1: The spin-structure function <j^ He (per nucleon) in the resonance region at Q 2 -values of 1.2, 
1.8, 2.4 and 3.0 (GeV/c) 2 . The error bars represent the total uncertainties with the inner part 



being statistical only. Also plotted are the DIS world data from E154 



2J and E99-117 



25( (note 



that these data are at dil erent Q 2 )- The curves were generated from the NLO parton distribution 
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to which target-mass corrections were applied. 



11 



0.005 



-0.005 



-0.01 



■ 


3 He 












BB 








GRSV 








AAC 


: o //' 






LSS 






o 


JLabE94-010 


- 

■ , , i , , , 


, i , 


■ 


This work 
, 1 , , 



-0.005 



-0.01 





neutron 




'/ 


.,§,..r.."'-'r' '* 




■ , , i , 


, , , 1 , , , 


, i , , 


i 


2 


3 




Q 2 (GeV/c) 2 





FIG. 2: T^ and Y\: test of spin duality on He (top) and neutron (bottom). Also plotted are 



the DIS parameterizations of Bliimlein and Bottcher 



17[ (grey band), GRSV 



181 ] (solid curve), 



AAC [19( (dashed curve) and LSS 
open circles are data from JLab E94-010 
uncertainty of each data set. 



201 ] (dotted curve) after applying target-mass corrections. The 
j|. The bands on the a>axis represent the systematic 
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FIG. 3: A[ in the resonance region. DIS data are from SLAC E142 



DESY experiment HERMES 



3J] and from JLab E99-117 



33(], E154 



241 ]. from 



251 ]. The error bars represent the total 



uncertainties with the inner part being statistical only. The curve represents a fit to the ^ He DIS 
data (see text). The arrows in the black frame point to the A(1232) peak position for each of our 
data sets 
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